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摘要 :【 目 的】 Wolbachia 是 一 种 广泛 存在 于 节肢 动物 中 的 胞 内 共生 细菌 ,影响 寄主 的 生物 学 特性 。 花 葡 马 
Frankliniella intonsa (Trybom) 是 重要 的 害虫 ,对 农作物 及 园林 植物 造成 危害 。 本 研究 旨 在 明确 Wolbachia THH Fy 
中 的 感染 情况 ,并 分 析 其 与 寄主 线粒体 DNA. 多 样 性 的 关系 。【 方 法 】 采 集中 国境 内 26 个 花 葡 马 自 然 种 群 ,运用 多 
位 点 序列 分 型 技术 (multilocus sequence typing, MLST) 对 其 体内 Wolbachia 感染 率 及 株 系 进行 分 析 ; 利 用 线粒体 COT 
分 子 标 记 研 究 花 蜀 马 的 遗传 分 化 及 遗传 多 样 性 ;通过 比较 感染 和 未 感染 Wolbachia 个 体 COI 数据 ,探究 Wolbachia 
多 样 性 与 寄主 线粒体 DNA 多 样 性 之 间 的 关系 。【 结 果 】 花 获 马 中 Wolbachia 的 感染 率 为 0% ~ 60% , 共 检 测 到 5 种 
Wolbachia 株 系 (Fintl ,wFint2 ,wFint3 ,wFint4 及 wFint5) , 均 属于 B 大 组 日 形成 一 个 单 系 群 。Wolbachia 感染 情况 与 
CHE FE tej Sh AP AE RE CC, GZ, TA 和 TY, N <5) 的 线粒体 DNA 多 样 性 相关 ,表现 为 不 感染 Wolbachia 的 种 群 中 线 粒 
体 DNA 单 倍 型 多 样 性 ( Hd) 与 核 蔡 酸 多 样 性 (Pi) 均 高 于 感染 Wolbachia 的 种 群 , 且 Wolbachia 感染 率 与 Hd 旦 显著 
负 相 关 (P<0.05)。AMOVA 分 析 表 明 花 欧 马 线粒体 DNA 遗传 分 化 与 Wolbachia 感染 情况 有 关 。[【 结论 ] Wolbachia 
可 能 在 侵 当 花 萄 马 种 群 后 出 现 遗 传 分 化 ; Wolbachia 感染 与 寄主 线粒体 DNA 多 样 性 有 关 。 

关键 词 : (EH; Wolbachia; 多 位 点 序列 分 析 ; 线粒体 DNA; 遗传 分 化 ; 遗传 多 样 性 
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Wolbachia infection and its relationship with mtDNA diversity in the 


























































































































flower thrips, Frankliniella intonsa ( Thysanoptera: Thripidae) 

LOU Heng, CHEN Ya-Ting, ZHANG Yan-Kai, YANG Xian-Ming, SUN Jing-Tao, HONG Xiao-Yue * 
( Department of Entomology, Nanjing Agricultural University, Nanjing 210095 , China) 

Abstract: [ Aim] Wolbachia is a bacterial endosymbiont that is found in a wide variety of arthropods, 
and affects the biological characteristics of host. The flower thrips, Frankliniella intonsa (Trybom) , is an 
important pest, and can cause damages to crops and garden plants. This study aims to examine the 
infection status of Wolbachia and its relationship with the mtDNA diversity of the host flower thrips. 
[Methods] The Wolbachia infection in 26 natural populations of F. intonsa was analyzed by using 
multilocus sequence typing ( MLST) in China. The genetic diversity and genetic differentiation of F. 
intonsa were investigated by using COI genetic markers. The COI data between the Wolbachia-infected 
and uninfected individuals were compared to explore the relationship between Wolbachia diversity and 
host mtDNA diversity. [Results] The prevalence of Wolbachia in the F. intonsa populations ranged from 
0% to 6096. Five Wolbachia strains (wFintl , wFint2, wFint3, wFint4 and wFint5) were detected. All 
the detected strains were assigned to supergroup B and formed a monophyletic group. The mitochondrial 
diversity of these populations (except CC, GZ, TA and TY, N <5) was associated with Wolbachia 
infection. The mtDNA haplotype diversity ( Hd) and nucleotide diversity ( Pi) of the uninfected group 
were higher than those in the infected group. Wolbachia infection was significantly negatively correlated 


with Hd (P «0.05). Furthermore, the analysis of molecular variance ( AMOVA ) revealed that the genetic 
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differentiation of F. intonsa mtDNA is related to the presence of Wolbachia infection. [ Conclusion ] 


Genetic differentiation of Wolbachia occurred after infecting host F. intonsa populations. Wolbachia 


infection is associated with host mtDNA haplotype diversity. 


Key words: Frankliniella intonsa; Wolbachia; multilocus sequence typing ( MLST) ; mtDNA; genetic 


differentiation; genetic diversity 


共生 菌 Wolbachia 4E BED Ics] vp oaa BR , 
其 主要 通过 母系 遗传 进行 传播 ( Hilgenboecker et 
al., 2008) 。 为 了 能 保证 在 寄主 种 群 中 稳定 遗传 ， 
Wolbachia 已 演化 出 多 种 调控 寄主 生殖 调控 的 机 制 ， 
包括 胞 质 不 亲 和 ( cytoplasmic incompatibility, CI) 、 
孤 肉 生殖 (parthenogenesis) , ASHE ( male-killing ) F0 WÈ 
性 化 (feminization ) 等 (Werren et al., 2008 ) ; 同时 ， 
Wolbachia 还 与 多 种 寄主 (如 秀丽 隐 杆 线虫 
Caenorhabditis elegans) 存在 共生 关系 (Darby et al., 
2012)。 一 些 研究 表明 , Wolbachia 还 参与 诱导 宿主 
免疫 系统 发 生 改 变 , 提高 对 病毒 及 其 他 病原 菌 的 抗 
HE, 如 Wolbachia 能 够 增强 黑 腹 果 蝇 Drosophila 
melanogaster 对 2 种 RNA 病毒 (Nora virus 和 Flock 
House virus ) 的 抗 性 (Teixeira et al., 2008) ,在 埃及 









































起 寄主 线粒体 DNA 多 样 性 降低 最 为 普遍 (Zhang et 
al., 2014) 。 基 于 Wolbachia 对 寄主 mtDNA 的 影响 ， 
在 应 用 mtDNA 进行 种 群 遗传 研究 时 要 充分 考虑 共 
生 菌 的 影响 因素 。 

44%] y Frankliniella intonsa ( Trybom ) Jgj ti] & 
科 ,是 一 种 通过 直接 吸食 汁液 而 破坏 多 种 观赏 性 植 
3] FE BY A Be Be AY AR Ex TE EL HR ( Hansen. et al., 
2003; Mound and Morris, 2007) ,也 是 多 种 植物 病毒 
的 传播 媒介 ( Ullman et al., 1997) 。 近 年 研究 表明 ， 
Wolbachia 侵 染 欧 马 科 昆 虫 , 且 能 够 调控 其 生殖 方式 
( Kumm and Moritz, 2008) ,如 Wolbachia 能 够 诱导 食 
A PES 5 EK eM Franklinothrips vespiformis RWE 
生殖 (Arakaki et al., 2001), 目前 Wolbachia EAE tij 
马 中 尚 无 相关 研究 。Wolbachia dé Jg& Ye dE tj 55? 


























伊 蚊 Aedes aegypti 中 , Wolbachia 能 够 抑制 登革热 病 
i fe Je ( Moreira et al.，2009 )。 近 年 来 , 关于 
Wolbachia 的 研究 引起 越 来 越 多 的 关注 。 

寄主 线粒体 DNA ( mitochondrial DNA, mtDNA) 
与 Wolbachia 一 起 营 母 系 遗 传 。 因 此 , Wolbachia 在 
寄主 中 的 广泛 分 布 与 传播 会 对 线粒体 DNA. 的 变异 
产生 影响 (Jiggins et al., 2003; Hurst and Jiggins, 
2005 )。 目 前 ,已 有 研究 表明 内 共生 菌 侵 染 与 寄主 
线粒体 DNA 关系 可 分 为 4 类 :(1) 内 共生 茵 驱动 的 
mtDNA 多 样 性 降低 ( Dean et al., 2003; Rasgon et al., 
2006; Raychoudhury et al.，2010 )。 例 如 ,感染 
Wolbachia HJ — 33 n d Tetranychus urticae 单 倍 型 及 
核 苷 酸 多 样 性 均 低 于 不 感染 个 体 (Yu et al., 2011) ; 
(2) 内 共生 戎 感染 不 引起 mtDNA 多 样 性 改变 。 例 
Till, Wolbachia 在 一 种 新 热带 区 甲虫 Chelymorpha 
alternans 中 的 感染 情况 与 寄主 单 倍 型 无 关 ( Keller et 
al., 2004) ; (3) 内 共生 菌 驱动 的 mtDNA 遗传 结构 的 
变异 。Baudry 等 (2003) 发 现在 绿 头 苍蝇 
Protocalliphora sialia 中 Wolbachia 株 系 与 特定 的 线 
粒 体 DNA 单 倍 型 相对 应 。Narita 等 (2006 ) 发 现 
Wolbachia 在 初始 不 感染 的 种 群 中 传播 将 会 导致 线 
粒 体 单 倍 型 频率 的 改变 ,影响 宽 边 黄粉 蝶 线 粒 体 遗 
传 结 构 ; (4) 内 共生 菌 驱 动 的 mtDNA 并 系 的 产生 
(Shoemaker et al., 2004), FL, Wolbachia 感染 引 




















花 葡 马 中 Wolbachia 能 否 影响 其 遗传 分 化 ?这 些 问 
题 对 于 研究 花 葡 马 种 群 动态 以 及 对 花 获 马 的 防治 等 
方面 具有 一 定价 值 。 本 研究 应 用 PCR 技术 检测 来 
自 中 国境 内 的 26 个 花 获 马 自然 种 群 Wolbachia 的 感 
染 情况 , 并 运用 多 位 点 序列 分 型 C multilocus 
sequence typing, MLST) 技术 明确 花 获 马 体内 
Wolbachia 的 株 系 类 型 ;通过 比较 分 析 花 萄 马 线粒体 
COI( cytochrome oxidase subunit I) 基因 与 Wolbachia 
PRA ITEM AA PRT T dE DUREE HP Wolbachia 
的 侵 染 历史 ;同时 ,比较 Wolbachia 的 感染 是 否 与 寄主 
线粒体 DNA 存在 相关 性 以 推断 Wolbachia 的 传播 模 
式 , 从 而 探究 Wolbachia 感染 对 花 获 马 线粒体 DNA 多 
样 性 的 有 影响。 研究 结 果 有 助 于 探究 Wolbachia 的 多 样 
性 及 其 与 宿主 遗传 结构 的 关系 ,同时 可 为 我 们 分 析 宿 
主 某 些 生物 学 特性 的 形成 机 制 \ 宿 主 遗 传 分 化 和 进化 
等 提供 重要 参考 信息 ( Rousset and Solignac, 1995) 。 





























1 材料 与 方法 


1.1 样本 收集 和 鉴定 

自 2009 年 6 H -2011 年 8 月 ,在 中 国境 内 共 采 
集 26 个 花 葡 马 种 群 ,共计 528 个 样本 ( 表 1, 图 1)。 
采样 时 ,在 花 葡 马 活动 的 植株 花 朱 下 放置 一 个 白色 
塑料 托盘 并 轻 拍 花 朱 , 或 将 花 激 单独 放 入 装 有 95% 
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R1 花 葡 马 样品 信息 及 Wolbachia 感染 情况 
Table 1 Information of Frankliniella intonsa samples and Wolbachia infection status 
种 群 代码 采集 地 点 地 理 坐 标 采集 时 间 样本 数 ”线粒体 单 倍 型 ” 感染 个 体 数 感染 率 (% ) Wolbachia 品系 
Population Collection Geographic Collection Sample mtDNA Infection Infection Wolbachia 
code locality coordinates date size haplotypes number rate strain 
北京 wFintl , wFint3 , 
BJ y 39.96?N, 116.33°E 2010.8 11 Hap28 , Hap45 4 36. 36 : 
Beijing wFint4 
BT AME 40.59eN, 109.84°E 2011.8 — 35 S pe 6 17.14 
Baotou, Inner Mongolia 33, 36, 45 
吉林 长 春 wFintl , wFint4 
CC a 43.88°N, 125.31°E 2011.8 4 Hap2, 6, 36 0 0 , d 
Changchun, Jilin : ape wFint5 
新 疆 昌 吉 Hap8 20 
CI 新 疆 昌 言 44.03°N, 87.39°E 2011.6 21 op a 1 4.76 wFint5 
Changji, Xinjiang 37, 38, 39 
云南 大 理 Hap30, 33, 40, wFintl , wFint3 , 
DL 25.60°N, 100.25°E 2011.7 45 1 42.22 
Dali, Yunnan 45, 46, 47 2 wFint4 , wFint5 
d E. 
GL J . 西 桂林 25.62°N, 110.67°E 2011.7 31 Hapl4, 16, 28 17 54. 84 wFintl , wFint4 
Guilin, Guangxi 
贵州 贵阳 Hapl2, 17, 19 
GY CANA ES 26.66°N, 106.82°E 2011.7 16 SR 4 25. 00 Fintl , wFint4 
Guiyang, Guizhou T 23, 24, 28 Wa ay 
广东 广州 n s ; 
GZ 23.16°N, 113.36°E 2011.7 2 Hap28 1 50.00 wFintl 
Guangzhou, Guangdong 
HA 江苏 淮安 33.66°N, 119.00°E 2009.6 29 Hapl6, 28, 45 14 48.28 wFintl , wFint4 
Huai’ an, Jiangsu 
黑龙 江 哈尔滨 Hap6, 7, 14, 32, 
HRB 45.74°N, 126.63°E 2011.8 32 4 12.50 Fint5 
E Harbin, Heilongjiang x 33, 36, 41, 45 TE 
1 iL Fintl , wFint3 
JJ . 江西 九江 29.53°N, 115.91°E 2011.7 10 Hap22, 28 6 60.00 eS 
Jiujiang, Jiangxi wFint4 
sey 
JQ ; 甘肃 酒泉 39.78°N, 98.51°E 2011.8 20 Hap33 0 0 
Jiuquan, Gansu 
云南 昆 昌 Hapl3, 16, 2 Fintl , wFint3 
KM 云南 屁 明 24.71°N, 102.72°E 2009.8 — 25 pe Mi 8 05.) ee eee 
Kunming, Yunnan 30, 42, 43, 45 wFint4 , wFint5 
—" 
LY 福建 龙岩 25.04°N, 116.43°E 2011.7 36 Hapl6, 19, 28 6 16.67 wFintl , wFint4 
Longyan, Fujian 
PET 
LZ HEM 36.05°N, 103.86°E 2011.8 24 Mp0, 33, 35, 0 0 
Lanzhou, Gansu 45, 48 
Ys ada SV Hapl6, 18, 1 
NP 福建 南 x 26.63°N, 118.22°E 2011.7 31 BELO TBs 12, 3 9.68 wFintl , wFint4 
Nanping, Fujian 26, 28, 45 
山东 青岛 E 5 : . 
QD : 36.32°N, 120.39°E 2009.6 6 Hap28 , 36 3 50.00 wFintl, wFint5 
Qingdao, Shandong 
3 Fa ik 
QHD 2 HIREA . 39.90°N, 119.54°E 2009.8 9 Hap28 , 45 4 44.44 wFintl , wFint4 
Qinhuangdao, Hebei 
EN 
SG 广东 韶关 24.68°N, 113.60°E 2011.7 10 Hap27 , 28 3 30. 00 wFintl 
Shaoguan, Guangdong 
河北 石家庄 Hapl6, 21, wFintl , wFint2 , 
SJZ 38.20°N, 114.36°E 2011.8 25 13 52.00 
J Shijiazhuang, Hebei d 28,45 wFint3 , wFint4 
i Tx Hap6, 14, 33 
SY 辽宁 沈阳 41.83°N, 123.57°E 2011.8 10 EE 2 20.00 wFint5 
Shenyang, Liaoning 36, 45 
山东 泰安 
TA 36.19°N, 117.12°E 2009.5 2 Hap6, 28 1 50.00 wFint4 
Tai’ an, Shandong 
TY 山西 太原 38.03°N, 112.65°E 2010.8 3 Hapl5, 36, 45 0 0 
Taiyuan, Shanxi 
江苏 泰州 Hapl, 3, 6, 16, 
TZ 32.37°N, 120.16°E 2010.7 46 18 39.13 Fintl , wFint4 
Taizhou, Jiangsu 25, 28, 45 2 wr 
青海 西宁 Hapl2, 31, 32 
XN . LEN x 36.60°N, 101.82°E 2011.8 13 P RE 4 30. 77 wFintl , wFint5 
Xining, Qinghai 33, 45 
吉林 延吉 H: 5,6, 10, 11 
YJ 吉林 延吉 42.89°N, 129.49°E 2011.8 32 pac duds 2 6.25 wFint5 


Yanji, Jilin 


14, 33, 34, 36, 44, 45 
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酒精 的 广 口 瓶 中 ,样本 间距 离 至 少 在 m 以 上 以 保 
证 样本 个 体 之 间 无 共同 的 亲本 。 

对 采集 的 样本 进行 鉴定 ,并 选择 花 蜀 马 进 行 下 
游 实验 。 花 葡 马 的 鉴定 主要 依据 成 虫 形 态 ,包括 头 部 
刚毛 的 位 置 . 大 小 和 数目 ,前 背 板 ,前 起 以 及 腹部 背 板 
工 ,颜色 特征 以 及 第 8 节 等 特征 (Funderburk , 2008 ) 。 
1.2 Wolbachia 检测 及 多 位 点 序列 分 型 

DNA 提取 :首先 向 已 编号 的 1.5 mL 离心 管内 
分 别 加 入 50 pL STE buffer ( 100 mmol/L NaCl, 10 
mmol/L Tris-HCl,1 mmol/L EDTA, pH 8.0) ,然后 分 
别 将 离心 管 置 于 冰 上 ; 挑 取 单 尖 花 葡 马 放 入 离心 管 
内 ,立即 用 玻璃 碾 梭 碾 磨 ; 碾 磨 后 将 离心 管 置 于 冰 上 
分 别 向 离心 管内 加 2 uL 和 蛋白酶 K(10 mg/mL); 
37°C 下 孵育 30 min; 孵 育 后 进行 953Y 变性 5 min, 以 
使 蛋白 酶 K 失 活 ;简单 离心 后 , -20Y 保存 或 立即 
作为 PCR 反应 的 模板 。 

Wolbachia 的 PCR 筛选 使 用 一 对 特异 性 引物 wsp- 
81F ( 5'-TGGTCCAATAACTGATCAAGAAAC-3') 和 
wsp-691R. (5'-AAAAATTAAACGCTACTCCA-3’) , 扩 增 
Wolbachia 表面 蛋白 ( Wolbachia surface protein) 基 
wsp 中 一 段 长 度 约 600 bp 的 序列 (Zhou et al., 1998) 。 

扩 增 条 件 为 :95Y 预 变性 3 min ;959C 30 s,52°C 
30 s,72C 1.5 min,35 个 循环 ;72% 延伸 7 min。 反 
应 中 包括 阴性 (确定 不 感染 Wolbachia 8f «3s 4E uj 
马 DNA) 和 阳性 对 照 (确定 感染 Wolbachia 的 单 头 花 
葡 马 DNA)。 通 过 1% 的 琼脂 糖 胶 进行 电泳 检测 , 凝 
胶 电 泳 确定 条 带 ,所 有 阳性 个 体 用 来 克隆 并 测序 
(上 海 美 吉 生 物 医药 科技 有 限 公司 ) 。 

切 胶 回收 目标 PCR 产物 ,通过 AxyPrep DNA AE 
胶 回 收 试剂 盒 (AxyCEN USA) 进行 产物 纯化 。 将 
纯化 产物 与 Trans-T1 载体 (TransGen Biotech ,北京 ) 
连接 ,再 转化 到 感受 态 大 肠 杆菌 Escherichia coli 中 ， 
将 阔 液 涂 布 于 含 50 mg/uL AIWA LB 培养 基 平 板 
上 37 倒置 恒温 培养 14 h, 观察 蓝 白 菌落 。 挑 取 白 
色 菌 落 置 于 已 加 入 1 mL 液体 LB 培养 基 ( 含 0.1% 
ATRAER ) 的 离心 管 中 ,37% 250 r/min 振荡 培养 
过 夜 ,然后 将 经 过 培养 后 的 菌 液 直 接 进行 PCR 检 
iM, £s PCR 检测 后 将 阳性 克隆 用 于 测序 。 

通过 分 析 wsp 序列 , 确认 花 蜀 马 感 染 一 种 
Wolbachia 株 系 ,随后 通过 5 个 用 于 Wolbachia MLST 
的 看 家 基因 (fisZ, coxA, fbpA, hcph 及 gatB) 进行 
PCR 扩 增 ,并 克隆 测序 ,具体 方法 见 Baldo 等 (2006 ) 
建立 的 分 型 方法 以 及 多 位 点 分 型 专用 数据 库 
(http://www. pubmlst. org/wolbachia/)。 看 家 基 






















































































序列 提交 至 Wolbachia MLST 数据 库 中 ,得 到 感染 花 
fij AY Wolbachia 序列 型 (sequence type, ST) 。 
1.3 453 COI 基因 扩 增 及 测序 

样本 通过 引物 C1-J-1751 和 C1-N-2329 扩 增 出 
COI 基因 , 约 571 bp 片段 (Simon et al. ,1994) , PCR 
反应 体系 为 25 kL, 包括 16.3 uL ddH,0,2.5 uL 
Dream Taq Buffer( 包含 2 mmol MgCl, ) ( Fermentas) , 
2.0 pL 2. 5 mmol dNTPs (TaKaRa) ,0.2 pL Dream 
Taq DNA 聚合 酶 (5 U/uL),1 uL DNA 以 及 上 下 游 
引物 各 0.5 uL(20 pmol/L) 。 扩 增 条 件 为 :95 民 预 变 
性 3 min;95°C 30 s,53% 30 s,72C 45 s,35 个 循环 ; 
TYCAERE 10 min。 扩 增 片段 通过 AxyPrep DNA XER 
回收 试剂 盒 (AxyCEN ,美国 ) 切 胶 回 收 并 进行 测序 。 
1.4 系统 发 育 分 析 

通过 Clustal X ( Thompson et al., 2002) 进行 序 
列 比 对 。 应 用 BioEdit( Hall, 1999) 进行 序列 的 编辑 
与 比 对 。 利 用 贝 叶 斯 推断 ( BE) 和 最 大 似 然 法 (ML) 
对 串联 的 MLST 基因 进行 系统 发 育 分 析 。 对 于 串联 
的 数据 ,通过 Jmodletest 选用 最 适 模 型 CTR +1+6 
模型 ,运用 该 模型 ,使 用 MrBayes 3. 1 进行 贝 叶 斯 分 
析 ( Ronquist and Huelsenbeck , 2003) ,分 析 从 初始 树 
开始 ,2 个 独立 的 运算 ,每 一 个 都 包括 4 个 链 , 运 算 
10 000 000 代 , 取样 频率 为 100, 并 且 最 初 的 
2 500 000 代 遗弃 。 为 剩余 树 计算 后 验 概率 。ML 树 
由 MEGA 5.0 软件 构建 (Tamura et al., 2011) 。 

COI 序列 使 用 DnaSP 5. 10 (Librado and Rozas, 
2009) 计算 的 单 倍 型 数 , 单 倍 型 多 样 性 (由 ) 和 核 苷 
酸 多 样 性 (Pi)。 在 Network 中 使 用 Median-joining 
( Bandelt et al., 1999) 对 COT 单 倍 型 种 内 的 系统 发 
育 进行 分 析 。 

1.5 ”种群 遗传 分 化 分 析 

通过 Arlequin v3.01 对 其 进行 分 子 变异 模型 分 
析 (AMOVA) ( Excoffier et al., 2005) 以 推断 线粒体 
DNA 基因 分 化 与 地 理 区 域 及 Wolbachia 感染 情况 的 
关系 。 以 长 江 为 界 , 将 花药 马 种 群 分 为 南部 种 群 
(TZ, JJ, NP, LY, SG, GZ, GL, GY, KM 和 DL) 和 
北部 种 群 (HA, QD, TA, TY, SJZ, BJ, QHD, SY, 
CC, YJ, HRB, BT, LZ, XN, JQ 和 CJ) ;根据 感染 
状态 ,将 每 个 地 理 种 群 的 花 蜀 马 拆 分 为 不 感染 亚 种 
群 和 感染 亚 种 群 ,并 将 种 群 分 为 感染 和 不 感染 两 大 
组 。 此 外 ,我 们 还 通过 SPSS 13.0 软件 对 单 倍 型 多 
样 性 ( 84) 与 22 个 种 群 ( 除 CC, GZ, TA MTY, N< 
5) Wolbachia 感染 率 进行 了 Pearson 相关 性 分 析 , 以 
探讨 Wolbachia 对 线粒体 遗传 多 样 性 的 影响 。 
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#2 2) Wolbachia MLST 等 位 基因 
2 结果 Table 2 MLST allelic profiles of Wolbachia detected from 


2.1 #23 BABE Wolbachia 感染 情况 

Wolbachia 检测 结果 表明 ,26 个 地 理 种 群 共 528 
个 个 体 中 , BR CC, JQ, LZ 和 TY 种 群 不 感染 
Wolbachia 外 ,其 余 22 个 地 理 种 群 的 花 葡 马 均 感 染 
Wolbachia。 检测 样本 中 共有 143 个 个 体感 染 
Wolbachia, 不 同 地 理 种 群 间 感 染 率 差异 较 大 ,在 
0% ~60% 不 等 ,平均 感染 率 为 27% ( 表 1) 。 
2.2 4E&j PRE Wolbachia 株 系 类 型 以 及 分 布 情况 

通过 对 Wolbachia wsp 序列 及 MLST 序列 的 分 
析 , 发 现在 26 个 花 葡 马 地 理 种 群 中 共 感 染 5 种 
Wolbachia 株 系 ,依据 命名 规则 命名 为 wFintl, 
wFint2, wFint3, wFint4 和 Fintg( 表 1)。 其 中 9 个 
地 理 种 群 (JL, JY, HA, LY, NP, QD, QHD, TZ 和 
XN) 发 现 双 感染 现象 。 通 过 对 单 感染 Wolbachia ( [XC 
有 1 条 特异 的 wsp 序列 ) 的 花 蓟 马 进行 多 位 点 序列 
分 型 技术 (MLST) 分 析 , 发 现在 感染 个 体 中 有 5 个 序 
列 型 (sequence types, STs) 与 wsp 序列 型 匹配 ( 表 
2)。 根 据 串 联 MLST 基因 系统 发 育 分 析 表 明 ,5 个 
STs 均 属 于 B 大 组 ,ST 间 遗 传 分 化 较 小 (图 1)。 

从 Wolbachia 株 系 的 分 布 来 看 , 22 个 感染 
Wolbachia 的 花药 马 种 群 中 ,17 个 地 理 种 群 ( BJ， 
BT, DL, GL, GY, GZ, HA, JJ, KM, LY, NP, QD, 




















Frankliniella intonsa 





E a E gatB coxA hcpÀ | ftsZ fopA 
wFintl 397 9 134 154 124 208 
wFint2 398 9 134 155 124 208 
wFint3 399 9 134 155 73 9 
wFint4 400 148 134 155 73 9 
wFintS 401 217* 201* 237* 124 9 





* MLST 数据 集 新 的 等 位 基因 Alleles new to the MLST database. 








QHD, SG, SJZ, TZ 和 XN) 检 测 到 wFintl ,感染 
率 约 为 50.3% ; 14 个 地 理 种 群 (BJ，BT，DL，GL， 
GY, HA, JJ, KM, LY, NP, QHD, SJZ, TZ ftl TA) 
中 检测 到 wFint4 ,感染 率 约 为 32.1% ; 仅 在 2 个 SJZ 
种 群 个 体 中 检测 到 wFint2 , 感染 率 极 低 (1. 39% ) ; 
在 5 个 种 群 中 (BJ，DL, JJ, KM 及 SJZ) 检测 到 
wFint3 ,其 感染 率 也 较 低 (2.2% ~10% ) ;在 9 个 种 
#E(BT, CJ, DL, HRB, KM, SY, QD, XN fl YJ) 中 
检测 到 wFint5 , 目地 理 分 布 较 广 ( 表 1, 图 2)。 

2.3 花 获 马 线粒体 COI 基因 单 倍 型 与 
Wolbachia 株 系 的 关系 

对 26 te tij = Ph FE (H 528 头 ) 进行 线粒体 
COL 基因 单 倍 型 分 析 ,我们 发 现 这 些 花 欧 马 种 群 中 
共有 48 个 单 倍 型 。 与 Wolbachia 感染 情况 相关 性 的 
分 析 发 现 , 感 染 Wolbachia 的 个 体 中 存在 15 种 单 倍 
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图 1 
Fig. 1 


Frankliniella intonsa ST397 

6|! Frankliniella intonsa ST398 
Frankliniella intonsa ST400 

Frankliniella intonsa ST399 


r— Acraea encedon ST3 





Armadillidium vulgare ST6 








0.99/83 Teleogryllus taiwanemma ST32 











Tribolium confusum ST30 











0.95/91 Drosophila simulans wNo ST16 








eek Ephestia kuehniella EkueB ST20 
Culex pipiens Cpip ST9 











0.99/95 





Protocalliphora sialia ST28 





0.99/96 — Encarsia formosa ST18 


Nasonia vitripennis ST26 


— Acraea eponina ST4 


us firmus ST21 


A Drosophila melanogaster ST1 


Drosophila simulans wRi ST17 





基于 MLST 看 家 基因 串联 序列 (2 079 bp) HARRY rf. Wolbachia 系统 发 育 分 析 


Phylogenetic of Wolbachia strains in Frankliniella intonsa based on the MLST concatenated data (2 079 bp) 
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图 2 Wolbachia 株 系 在 中 


@ wFintl 
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© 未 感染 Uninfected 


国花 蜀 马 种 群 中 的 分 布 











Fig. 2 Distribution of Wolbachia variants in the examined populations of Frankliniella intonsa collected across China 
地 图 来 自 国家 测绘 局 官方 网 站 1 亿 示 意 地 图 版 3。Map from SBSM 100 million cartogram 5. 





型 ,不 感染 个 体 中 存在 43 种 单 倍 型 。 相 同 线粒体 
DNA 单 倍 型 可 以 携带 不 同 Wolbachia 株 系 ,如 单 倍 
型 H28,H45 和 H30 分 别 与 4, 3 和 2 个 Wolbachia 
株 系 相对 应 。 反 之 亦 然 ,相同 的 Wolbachia 株 系 也 与 
不 同 COL 单 倍 型 相对 应 ,如 wFintl 在 5 个 单 倍 型 出 
现 (H25,H28,H30,H41 及 H45) , wFintS 在 6 个 单 
倍 型 中 出 现 (H32, H33, H36, H37, H40 及 H45), 
wFint2 仅 与 一 个 单 倍 型 对 应 , 且 只 存在 于 一 个 种 群 
(SJZ) (图 2 和 3)。 
2.4 Wolbachia 感染 与 花 葡 马 线粒体 DNA 多 样 性 
有 关 

研究 表明 ,感染 Wolbachia 的 花 萄 马 种 群 单 倍 型 
多 样 性 (Hd) ZERE FEE (Pi) 显著 低 于 未 感染 
种 群 ( 感 染 种群 : Hd =0.566, Pi =0.00168; 不 感染 
种 群 : Hd = 0. 846, Pi = 0. 00435) ( 表 4)。 同 时 ， 
Pearson 相关 分 析 表 明 花 葡 马 单 倍 型 多 样 性 ( 4) 与 
Wolbachia 感染 率 呈 显著 负 相 关 性 (R= -0.556”， 
P<0.05, 图 4)。 

种 群 遗传 分 化 分 析 ( AMOVA) 结果 显示 无 论 以 
哪 种 方式 分 组 ,组 之 间 的 分 化 都 达到 了 极 显著 的 水 
平 (P<0.0001)( 表 5)。 数 据 集 1(dataset 1) 显示 
不 同 地 理 区 域 单 倍 型 之 间 存 在 显著 分 化 ,表明 在 这 
些 种 群 由 于 存在 较 强 的 地 理 隔离 ,而 存在 较 小 的 基 
因 流 。 数 据 集 2( dataset 2) 显示 单 倍 型 和 感染 状态 
ZAMIR (Fy =0.39816,P <0.0001) ,这 一 结 
表明 花 葡 马 较 高 的 遗传 分 化 发 生 在 不 同感 染 类 型 的 
个 体 间 ,同时 , 花 萄 马 线粒体 单 倍 型 并 不 时 随机 分 
布 ,而 是 与 其 Wolbachia 感染 类 型 相关 。 在 两 种 分 组 








方式 中 ,按照 感染 状态 分 组 方式 下 的 组 之 间 的 变异 
所 占 的 比例 比 按照 地 理 区 域 分 组 方式 下 的 组 之 间 的 
比例 高 , 达到 9. 2496 , 也 就 是 说 花 获 马 个 体 间 
Wolbachia 感染 与 否 造成 的 遗传 分 化 比 地 理 区 域 造 
成 的 遗传 分 化 更 为 显著 ,同时 证 明了 线粒体 DNA 的 
遗传 分 化 与 Wolbachia 存在 一 定 的 联系 ( 表 5) 。 


3 讨论 


本 研究 检测 26 MEH HAPEE Wolbachia 感染 率 
为 0% ~ 60% ,说 明 Wolbachia FEW IL UY deg Hy, PA 
I FPA Ee A 4E 8j SS BRE Wolbachia 感染 率 都 为 
100% ,不 同 地 理 种群 Wolbachia 感染 率 仍 有 较 大 差 
5t. HBT, Wolbachia 在 花 葡 马 中 的 传播 率 对 寄主 适 
合 度 的 影响 以 及 CI 等 生殖 方面 的 影响 都 需要 进 一 
步 的 研究 。 

通过 分 析 花 蜀 马 种 群 mtDNA 序列 及 感染 
Wolbachia 的 株 系 ,可 揭示 二 者 之 间 的 内 在 联系 。 基 
F wsp 序列 及 MLST 分 析 发 现 5 种 Wolbachia IRA, 
且 株 系 间 wsp 序列 变异 较 小 。 这 一 结果 表明 ， 
Wolbachia 株 系 可 能 共同 来 自 于 历史 上 某 一 次 侵 染 
的 一 种 Wolbachia 株 系 ，Wolbachia 可 能 初始 侵 染 后 
产生 遗传 分 化 。 研 究 还 发 现 感染 及 未 感染 
Wolbachia 个 体 间 存在 共享 单 倍 型 ,造成 这 种 现象 的 
原因 可 能 是 当 Wolbachia 垂直 传播 率 不 高 时 ， 
Wolbachia 在 世代 传递 中 丢失 ,但 是 由 于 mtDNA 和 
Wolbachia 同时 营 母 系 遗 传 (Hinrich et al., 2002; 
Charlat et al., 2009) ,即使 在 Wolbachia 感染 丢失 的 
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#3 Wolbachia 感染 情况 以 及 不 同 种 群 中 单 倍 型 分 布 


Table 3 Infection rate and haplotype diversity of Wolbachia amongst different populations of Frankliniella intonsa 





图 














A HOH were 
种 群 代码 P kem 线粒体 DNA 单 倍 型 eee E d 线粒体 DNA 音信 型 
Population code gei mtDNA haplotype un mtDNA haplotype 
individuals individuals 
BJ 4 Hap28 7 Hap28, 45 
BT 6 Hap6, Hap36, Hap45 29 Hap2, 6, 14, 33, 36, 45 
CC 0 4 Hap2, 6, 36 
CJ 1 Hap37 20 Hap8, 9, 20, 37, 38, 39 
DL 19 Hap30, 33, 40, 45, 46 26 Hap30, 45, 46, 47 
GL 17 Hap28 14 Hapl6, 28 
GY 4 Hap23 28 12 Hap17, 19, 23, 24, 28 
GZ 1 Hap28 1 Hap28 
HA 14 Hap28 15 Hapl6, 28, 45 
HRB 4 Hap36, Hap45 28 Hap6, 7, 14, 32, 33, 36, 41, 45 
JJ 6 Hap22, 28 4 Hap28 
JQ 0 20 Hap33 
KM 8 Hap30, 45 17 Hap13, 16, 29, 30, 42, 43, 45 
LY 6 Hap28 30 Hapl6, 19, 28 
LZ 0 24 Hap6, 33, 35, 45, 48 
NP 3 Hap28 28 Hapl6, 18, 19, 26, 28, 45 
QD 3 Hap28 , 36 3 Hap28 
QHD 4 Hap28 , 45 5 Hap28 
SG 3 Hap28 7 Hap27 , 28 
SJZ 13 Hap21 , 28 12 Hapl6, 28, 45 
SY 2 Hap36 8 Hap6, 14, 33, 36, 45 
TA 1 Hap28 1 Hap6 
TY 0 3 Hap15, 36, 45 
TZ 18 Hap25, 28, 45 28 Hapl, 3, 6, 16, 25, 28, 45 
XN 4 Hap32, 45 9 Hap12, 31, 32, 33, 45 
YJ 2 Hap36, 41 30 Hap4, 5,6, 10, 11, 14, 33, 34, 36, 44, 45 
总 计 Total 143 15 385 43 
Q wFintl 
Q wFin2 
@ wFin3 
Q wFin4 
® wFint5 
© 未 感染 Uninfected 


3 FEB 5 REIS DNA 单 倍 型 中 接 网 络 分 析 及 Wolbachia 感染 情况 


Fig. 3 Median-joining networks of Frankliniella intonsa mtDNA ( COI) haplotypes associated with Wolbachia infection status 
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X4 花 葡 马 线粒体 DNA 单 倍 型 多 样 性 及 核 彰 酸 多 样 性 


Table 4 mtDNA and nucleotide diversity in Frankliniella intonsa 














种 群 代码 感染 情况 单 信 型 数量 AIRE BFE REE 
Population code Infection status Nh Hd Pi 
BI wFintl , wFint3, wFint4 2 0. 182(0. 144) 0. 00032 (0. 00025) 
BT wFintl , wFint , wFintS 6 0. 765(0.045) 0. 00415 (0. 00052) 
CC 3 0. 833 (0. 222) 0. 00668 (0. 00192 ) 
CJ wFint5 6 0.638(0.098) 0.00425 (0. 00070 ) 
DL wFintl , wFint3 , wFint4, wFint5 6 0.541(0.072) 0.00174(0. 00045) 
GL wFintl , wFint4 2 0.125(0.077) 0.00087 (0. 00054 ) 
GY wFintl , wFint4 5 0. 700(0. 080) 0. 00530 (0. 00050 ) 
GZ wFintl 1 0.000(0. 000) 0. 00000 (0. 00000 ) 
HA wFintl , wFint4 3 0. 135(0. 085 ) 0.00059 (0. 00044 ) 
HRB wFint5 8 0. 839(0.032) 0.00379 (0. 00050 ) 
JJ wFintl , wFint3, wFint4 2 0.200(0. 154) 0. 00035 (0. 00027 ) 
JQ 1 0. 000(0. 000) 0. 00000 (0. 00000 ) 
KM wFintl , wFint3, wFint4, wFint5 7 0. 580(0. 108) 0. 00182(0. 00054) 
LY wFintl , wFint4 3 0. 560(0. 068) 0.00408 (0. 00038 ) 
LZ 5 0.667 (0.063) 0.00225 (0. 00065 ) 
NP wFintl , wFint4 6 0.619(0.080) 0.00424 (0. 00050 ) 
QD wFintl , wFint5 2 0.333(0.215) 0.00116(0. 00075 ) 
QHD wFintl , wFint4 2 0.222(0. 166) 0.00039 (0. 00029 ) 
SG wFintl 2 0.200(0. 154) 0.00035 (0. 00027 ) 
SJZ wFintl , wFint2, wFint3, wFint4 4 0.230(0. 110) 0. 00096 (0. 00055) 
SY wFint5 5 0. 800 (0. 100) 0. 00329 (0. 00103 ) 
TA wFint4 2 1.000(0. 500) 0. 00871 (0. 00436 ) 
TY 3 1.000(0. 272) 0.00581 (0. 00226 ) 
TZ wFintl , wFint4 7 0. 386(0. 090 ) 0. 00232 (0. 00065 ) 
XN wFintl , wFintS 5 0.756(0. 097) 0. 00241 (0. 00070 ) 
YJ wFint5 13 0.821 (0.061) 0. 00479 ( 0. 00059 ) 
IN 感染 个 体 Infected individuals 15 0. 566(0.046) 0.00168(0.00135) 
UN 不 感染 个 体 Uninfected individuals 43 0.846(0.011) 0. 00435 (0. 00242 ) 
括号 中 数值 示 标 准 差 。Values in parentheses are SD. 
1.000 
R--0.556 P«0.05 
0.800 BA o 
o 
= o 
Š 0.600 Š a 
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感染 率 Infection rate (96) 
IK 4 Pearson 相关 系数 分 析 单 倍 型 多 样 性 与 Wolbachia 感染 率 之 间 的 关系 











Fig. 4 Pearson correlations between haplotype diversity (Hd) and the infection rate of Wolbachia in Frankliniella intonsa 
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表 5 AMOVA 分 析 花 葡 马 遗传 多 样 性 
Table 5 AMOVA analysis of the genetic variation in Frankliniella intonsa 
"deo TE 变异 来 源 由 度 平方 和 变异 组 分 占 总 变异 比例 (% ) 固定 系数 P 值 
Source of Sum of Variance Percentage of Fixation 
Dataset Subsets MN df rh : P-value 
variation squares component variation index 
数据 集 1( 所 有 个 体 ,地 理 区 域 分 组 ) 
Dataset 1 (all, divided into geographical regions) 
组 间 恋 屁 
分 组 癌变 民 1 27.562 ”0.06303Va 5.46 
Among groups 
乡 "Ep c 
组 内 种 群 癌变 异 24 182.048 0.34862Vb 30. 18 
South, Among populations within groups 
1 North 种 群 内 变异 0.35640 <0.0001 
i HM 503 — 373.915 0.74337Vc 64.36 
Within populations 
P4 zl i 
d 528 583.526 1.15502 
Total 
数据 集 2( 所 有 个 体 ,感染 情况 分 组 ) 
Dataset 2 (all, divided according to infection status ) 
PET 
HER 1 28.258 0. 10905Va 9.24 
Among groups 
2] $ alar c 
组 内 种 群 癌变 异 46 212.332 0.36083Vb 30. 58 
Infected Among populations within groups 0. 39816 0.0001 
Uninfected 种 群 内 变异 ` S. 
AP . 483 343.043 .0.71023Vc 60.18 
Within populations 
P4 Ap i 
SER 530 583.633 1.18011 
Total 


情况 下 ,与 其 相关 联 的 mtDNA 单 倍 型 仍 能 够 随 着 母 
体 细胞 遗传 给 不 感染 个 体 ,通过 这 个 过 程 的 不 断 累 


W ,不 感染 个 体 获得 与 感染 个 体 一 致 的 单 倍 型 。 另 


外 ,我 们 还 通过 检测 花 荀 马 种 群 mtDNA 序列 来 分 析 
Wolbachia 感染 与 寄主 个 体 mtDNA 多 样 性 及 进化 有 
Ko 不 感染 种 群 Hd 和 Pi 比 感染 种 群 高 , 同时 ， 
Pearson 相关 分 析 表 明 Hd Bt Wolbachia 感染 率 的 增 
加 而 下 降 。 基 于 以 上 两 点 ,推断 Wolbachia 的 感染 可 
能 导致 线粒体 DNA 多 样 性 下 降 。 同 时 ,通过 分 子 变 
异 模 型 分 析 , 进 一 步 验 证 了 花 萄 马 种 群 遗传 分 化 与 














进行 遗传 分 化 分 析 时 ,除了 考虑 地 理 因 素 以 外 ,还 须 
考虑 寄主 Wolbachia 感染 情况 。 我 们 推测 Wolbachia 
感染 引起 遗传 分 化 可 能 是 由 于 其 对 寄主 生殖 的 影响 。 
Wolbachia 通过 对 寄主 产生 生殖 调控 ,导致 寄主 生殖 
隔离 ,影响 寄主 种 内 基因 交流 ,从 而 导致 寄主 线粒体 
DNA 多 样 性 降低 (Dean et al., 2003; Rasgon et al., 
2006; Raychoudhury et al., 2010; Smith et al., 2012) 

通过 研究 分 析 花 萄 马 中 Wolbachia 株 系 的 系统 
发 育 关系 及 侵 染 历史 ,揭示 了 Wolbachia f AES = 
种 群 并 且 出 现 遗 传 分 化 。 同 时 , Wolbachia 感染 可 能 



































Wolbachia 的 感染 类 型 相关 。 与 我 们 结果 类 似 , 库 蚊 
中 Wolbachia 感染 降低 寄主 的 线粒体 多 样 性 (Rasgon 
et al., 2006), Raychoudhury & (2010) 在 对 金 小 蜂 
的 研究 中 也 发 现 ,Wolbachia 感染 对 寄主 线粒体 存在 
一 种 选择 性 扫除 效应 (selective sweep) ,从 而 降低 了 
宿主 的 线粒体 多 样 性 。Wolbachia 感染 的 一 种 果 蝇 
Drosophila recens 的 线粒体 多 样 性 比 不 感染 
Wolbachia 姐妹 种 Drosophila subquinaria 的 线粒体 多 
样 性 低 (Shoemaker et al., 2004)。 同 时 ,发 现 截 形 叶 
螨 的 Wolbachia 感染 个 体 比 不 感染 个 体 线粒体 DNA 
单 倍 型 多 样 性 ( Hd) RECTE BR AE FEE CPi) 低 ( Zhang 
et al., 2013) 。 

研究 还 表明 , Wolbachia 感染 较 地 理 区 域 对 花 敬 
马 种 群 遗 传 分 化 的 影响 更 为 显著 。 因 此 ,在 对 寄主 














导致 线粒体 DNA 多 样 性 降低 。 我 们 之 后 还 将 从 花 
葡 马 中 维持 Wolbachia 的 感染 机 制 与 适应 性 ,垂直 传 
播 率 以 及 胞 质 不 亲 和 强 度 ,其 与 线粒体 DNA. 互 作 机 
制 等 其 他 方面 进行 进一步 研究 ,从 而 揭示 Wolbachia 
与 线粒体 DNA 之 间 的 关系 以 及 Wolbachia —; 468i = 
的 互 作 关系 。 
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